NASA Technical Memorandum 105213 


// 7 - 

- -3 

if 

p 

Differential Continuum Damage Mechanics Models 
for Creep and Fatigue of Unidirectional 
Metal Matrix Composites 

S.M. Arnold 

Lewis Research Center 

Cleveland, Ohio 

and 

S. Kruch 

Office National D ’ Etudes et de Recherches Aerospatiales 
Chatillon, France 


November 1991 






DIFFERENTIAL CONTINUUM DAMAGE MECHANICS MODELS FOR 

CREEP AND FATIGUE 

OF UNIDIRECTIONAL METAL MATRIX COMPOSITES 

By 

S. M. Arnold 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland Ohio, 44135. 

and 

S. Kruch 

Office National d’ Etudes et de Recherches Aerospatiales 
92322 Chatillon, France. 


ABSTRACT 

Three multiaxial isothermal continuum damage mechanics model for creep, 
fatigue, and creep/fatigue interaction of a unidirectional metal matrix composite 
volume element are presented, only one of which will be discussed in depth. Each 
model is phenomenological and stress based, with varying degrees of complexity to 
accurately predict the initiation and propagation of intergranular and transgranular 
defects over a wide range of loading conditions. The development of these models 
are founded on the definition of an initially transversely isotropic fatigue limit 
surface, static fracture surface, normalized stress amplitude function and isochronous 
creep damage failure surface, from which both fatigue and creep damage evolutionary 
laws can be obtained. The anisotropy of each model is defined through physically 
meaningful invariants reflecting the local stress and material orientation. All three 
transversely isotropic models have been shown, when taken to their isotropic limit, to 
directly simplify to previously developed and validated creep and fatigue continuum 
damage theories. 

Results of a nondimensional parametric study illustrate i) the flexibility of the 
present formulation when attempting to characterize a large class of composite 
materials and ii) its ability to predict anticipated qualitative trends in the fatigue 
behavior of unidirectional metal matrix composites. Additionally, the potential for 
the inclusion of various micromechanical effects (e.g. fiber/matrix bond strength, 
fiber volume fraction, etc.), into the phenomenological anisotropic parameters are 
noted, as well as a detailed discussion regarding the necessary exploratory and 
characterization experiments needed to utilize the featured damage theories. 


1.0 INTRODUCTION 

In recent years, due to advances in scanning electron microscopy, acoustic 
emission techniques as well as other nondestructive testing methods, a consensus has 
emerged that the nonlinear response of solids and ultimately their mechanical strength 
are dependent not only on the basic structure of the material but also on the type, 
distribution, size and orientation of the defects in its structure. As a result a 
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relatively new branch of continuum mechanics, known as Continuum Damage 
Mechanics (CDM), has emerged. Like Continuum Mechanics, which allows one to 
describe (over an appropriate representative volume element) the heterogeneous 
microprocesses involved during the straining of materials and structures at the 
macroscale; CDM allows one to describe the material’s progressive deterioration 
(damage) from the virgin state (no damage) to the final state, corresponding 
generally to macro-crack initiation (or the "breaking up" of the representative volume 
element). Therefore, CDM allows one to predict the life limiting macroscopic 
properties such as rupture strength, fatigue life, or creep rupture lifetime of the 
material. 

The tracking, or description, of the evolution of damage is accomplished through 
the introduction of special thermodynamic (internal) field variables representing in an 
appropriate statistical sense, the distribution and density of defects locally. 

Numerous damage theories, both micromechanical [1—6] 1 and phenomenological [7-11], 
for example, have been proposed and discussed in the literature. The reason behind 
such diversity in the mathematical nature of the damage variable(s) (e.g., scalar 
[7,12,13], vectors [14,15] and tensors [12,16,17]) and thus the damage theories, stems 
from the difficulty associated with directly measuring "damage" macroscopically and 
the degree of approximation with which the internal variables describe the salient 
aspects of the macroscopic effects of the micro-defect kinematics. For a number of 
excellent review articles and books on this subject the reader is referred to references 
[7,12,18-24]. 

In this paper a CDM isothermal transversely isotropic creep, fatigue, and 
creep/fatigue model, with a scalar damage internal variable, oriented toward 
unidirectional metallic composites, will be described. These models are extensions of 
previously developed and validated models, developed at ONERA (Office Nationale 
d’Etudes et de Recherches Aerospatiales), [10,11,18,19,25,26,27,28] for isotropic 
monolithic metals. As a result, the paper begins with a brief review of CDM in 
general, and in particular, the work undertaken by Chaboche and his colleagues at 
ONERA, so as to provide a foundation to describe the extension of the ONERA 
models to transversely isotropic materials (e.g., metal matrix composites). 

Subsequent to the multiaxial statement of the transversely isotropic CDM models, a 
nondimensional uniaxial parametric study is undertaken to illustrate the various 
features and flexibilities of the proposed extensions. The paper then concludes with a 


‘Square brackets, i.e. [ ], will indicate suggested references for further reading. 
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discussion regarding exploratory and characterization experiments required to utilize 
and validate the proposed forms. 

2.0 MEASURES AND DEFINITIONS OF DAMAGE 

Material (lattice) defects can be roughly classified with respect to their 
geometry into i) point defects; for example vacancies, interstitial and impurity 
atoms, ii) line defects (dislocations), iii) plane defects (slip planes and cracks), and 
iv) volume defects (cavities and inclusions). Numerous ways are available to define 
the internal variables associated with these defects and damage processes. Each 
definition, however, must correspond to some method of measurement and implies a 
given approach, either a micromechanical or phenomenological one. Here, four typical 
approaches will be briefly stated. 

Measurements at the microstructural scale (i.e., density of micro-cracks in 
fatigue or the volume fraction of cavities in ductile damage) lead to microscopic 
models that can be integrated over the macroscopic volume element to obtain 
the properties of the damaged volume element. However, difficulty is 
encountered when incorporating these results into life prediction methods. 

Measurements of macro physical parameters (such as density, acoustic emission 
and resistivity) may lead to some consistent macroscopic definitions of damage 
parameters but require the definition of a global model to convert them into 
properties which characterize mechanical resistance. 

Measures of remaining life are used to account for the cumulative aspects of 
damage processes in the life prediction methods but do not directly lead to a 
damage constitutive law. They do however, provide insight into interesting 
properties of damage, for example in fatigue they indicate that the damage 
evolution equations must be non-separable relative to the damage and loading 
variables [29]. 

Measurements of variations in the macro mechanical behavior (i.e., modification 
of elastic, plastic or viscoplastic properties) are the most appropriate from a 
phenomenological point of view. These measurements are easier to interpret in 
terms of damage variables using the concept of effective stress [12,28] with an 
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equivalence in strain, Fig. 1. This is the approach that will be followed 
throughout the remainder of this report. 

The effective stress (cr) concept, with an equivalence in strain as put forth by 
Lemaitre and Chaboche [12] states that a damaged volume of material under the 
applied stress (u) shows the same strain response (whether the behavior is elastic, 

plastic, or viscoplastic) as the undamaged one submitted to the effective stress (a). 
In mathematical form 


a 



a 

T=U 


( 1 ) 


where A = A-A^ and the damage D represents the loss of effective area (A^/A) 
taking into account decohesions and local stress concentrations. This definition has 
been supported by the results of homogenization techniques [30] and damage measures 
utilizing this concept have been demonstrated under various situations, for example, 
in the case of ductile rupture [31], creep damage [32], and fatigue damage for both 
monolithic and composite materials [18,33,34]. 

3.0 CDM MODELS FOR INITIALLY ISOTROPIC MATERIALS 


Different damage variables are associated with different damage processes, such 
as creep, fatigue, ductile and brittle damage, and have been discussed in the 
literature extensively as noted above. Here we will confine ourselves to the modeling 
of creep, fatigue and creep /fatigue damage using a scalar damage measure in 
monolithic metals and then extend this discussion, in subsequent sections, to metal 
matrix composites by allowing the anisotropic evolution of this measure. Figure 2 
clearly illustrates the damage mechanisms for the three cases of interest, i.e., creep, 
fatigue and creep/fatigue interaction. 

3.1 Creep Damage 

In Fig. 2a the two primary damage mechanisms [35,36], nucleation and growth 
of intergranular defects (typically by a diffusion process) are illustrated schematically 
when the material is subjected to a pure creep loading condition. Recently, Hayhurst 
[37] has shown that it is possible to make some connection between the equations 
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obtained from a materials science approach and the more macroscopic ones developed 
under the framework of CDM. 

CDM was first developed for the case of creep damage by Kachanov [38] and, 
Rabotnov [23], since then the concept of effective stress has been shown to predict 
tertiary creep curves as well as changes in creep ductility [39]. Leckie and Hayhurst 
[8,13,40,41] have generalized the classical [23,38] uniaxial creep damage equations to 
multiaxial stress conditions by describing isodamage surfaces (or isochronous surfaces) 

defined in terms of three stress invariants i) the octahedral shear stress which 

is related to the effects of shear, ii) the hydrostatic stress J^{cr), which greatly 
affects the growth of cavities; and iii) the maximum principal stress Jq(ct) = ^ max > 
which opens the micro-cracks and causes them to grow. The equivalent stress is 
then defined through a linear combination of these invariants : 


X(a) = a c J Q (a) + 0 c ^(ff) + (\-a~PJ J 2 (a) 


( 2 ) 


here and are coefficients dependent upon the material. The form of the 
resulting creep damage model is given in Table I. Note that here the exponent k is 
taken to be a constant, thus implying linear accumulation of creep damage even 
though the evolution of damage may be nonlinear. 

Furthermore, the results of creep-like cyclic tests (e.g., high frequency-low 
amplitude or high amplitude-low frequency) have shown the need to include some 
factor relating the response time of the material versus that of the rate of loading. 
For example, when a material is subjected to a high frequency loading with a mean 
stress equal to zero, the material damage induced by creep is negligible, thereby 
leading to the case of "pure" fatigue damage. However, when the material is 
subjected to a high frequency and a low or medium amplitude with a high mean 
stress, the damage induced by creep is dominant. Therefore a new material 


parameter (r) and a "delayed" stress (cr?.) were introduced into the formulation, so 
as to account for the lag time between load application and material response [26]. 
As a result, given a high frequency load history, the delayed stress is equal to the 
mean stress while for a low frequency load history the delayed stress is the real (or 
applied) stress. 
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3.2 Fatigue Damage 


In Fig. 2b, the fatigue crack initiation and growth process is illustrated 
schematically, with the damage (D) being associated , macroscopically, with the 
initiation and propagation of transgranular defects (e.g., slip bands and micro-cracks). 
For example the micro-crack measurements made by Cailletaud and Levaillant [42], 
Hua and Socie [43], and Socie et. al. [44], have shown the possible equivalence 
between the definition of D by the effective stress concept, the definition of D in 
terms of the remaining life concept, and the quantification of physical damage in 
terms of micro-cracking. A number of fatigue damage models have been proposed in 
the past [12,18,21]. A general form, depending upon one’s choice for a (the exponent 
which may also be a function of stress), as indicated by Chaudonneret and Chaboche 
[45], that leads to rules considered earlier by various authors [10,25,26,46,47,48] is 
given by: 


dD = D^ a m’^ 


o '—a 
m 1 


{ mitt 


dN 


( 3 ) 


where a and ~d are, respectively, the maximum and mean stress within a cycle N. 

In the case of fatigue, several aspects need to be considered when developing a 
phenomenological model: 


a) the existence of micro-initiation and micro-propagation stages, 

b) the existence of a fatigue or endurance limit and its marked decrease 
after prior damage, 

c) the existence of a static limit (ultimate stress level), 

d) the effect of mean-stress on the fatigue limit or the S-N curves 
themselves, and 

e) the typical nonlinear-cumulative effects for two level tests or 
block-program loading conditions. 

Three CDM fatigue models, developed at ONERA, which to varying degrees 
accurately predict the above observations, are given in multiaxial form in Table II 
and uniaxial form in Table III. Note that these models are obtained from 
equation (3) (predominantly developed from a remaining life concept) and 
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incorporated into a CDM model utilizing the effective stress concept, by a convenient 

* 

change of variable (i.e. , D ), that is : 


D = 1 - (1 - D*) 1+/? (4) 

The ability of equation (3), and therefore those fatigue models in Table II and 
III, to represent the required nonlinear accumulation during a two level or 

block-program loading test is directly linked to the dependency of a on a m and a ; 
thus rendering the equations non-separable. It is this nonseparability and not the 
nonlinearity of the damage evolution equation that permits the modeling of a 
nonlinear accumulation of damage [29]. Model I, was proposed in 1974 [49] and has 
since been exercised for a number of different applications and materials [50-52], 
Results indicate good predictive capabilities under varying fatigue cycles, even though 
a single scalar state variable D represents both micro-initiation and 
micro-propagation periods. However, due to this lack of separation, a deficiency in 
Model I was observed when analyzing fatigue loading histories with prior creep or 
creep like (e.g., high cycle, low amplitude fatigue cycles) loading histories. 

Furthermore microstructural observations have indicated that fatigue damage can be 
divided into two periods, [53,54] the initiation period (corresponding to the creation 
of micro-cracks) and the propagation period (which is characterized by the inward 
growth in the material of a few micro-cracks). 

Therefore in Model II a separation between the micro-crack initiation and 
propagation stages was introduced. Note that in Model II, a is assumed to be 
constant, thereby producing only linear damage accumulation (Miner’s rule) during 
the propagation stage. This does not, however, imply that the damage D (i.e., D = 

initiation + propagation) is restricted t0 onl y linear accumulation, as the damage 
accumulation can be shifted by an increase or decrease in the number of cycles for 
crack initiation. As a result model II can be shown to give the Manson Double 
Linear Damage Rule (OLDR). Additionally, in model II, the initiation of 
microcracks are influenced by the frequency (v) of the applied loading due to the fact 
that differences in the experimental results were observed when varying the frequency 
of loading, particularly when at low frequency (e.g., 0 < v < 10 Hz) levels. 

Model III has been recently suggested as a compromise between Models I and 
II, in that it has the simplicity of Model I (no separation between initiation and 
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propagation is assumed) yet gives rise to a non evanescent fatigue limit by the 
introduction of a term representing the fatigue limit of propagation (ct^). 

3.3 Creep-Fatigue Interaction 

In Fig. 2c the interaction of intergranular and transgranular defects is depicted 
schematically. It is surmised, within the pertinent domain where creep-fatigue 
interaction is applicable, that the presence of cavities allow for easier crack propaga- 
tion and that the increase in stress intensity at a crack tip causes an increase in the 
nucleation and coalescence of voids. 

Analytically this interaction is represented using the effective stress concept, by 
assuming that the mechanical effects of creep and fatigue damage can be directly 
added, i.e., 


dD = dD c + dDp 

= f c (a,T,D,...) dt + f F (a m ,a,T,D,...) dN (5) 

where f and fp represent the characteristic functions of creep and fatigue respectively 
as described above and D=D c +Dp. Thus the two functions (f c and fp) can be 
determined independently from pure tensile creep tests and pure high frequency 
fatigue tests. The conditions under which accumulation of both macroscopic effects 
would take place (e.g., low frequency or loadings with hold times) are then predicted 
by integrating numerically the above equations. This approach has been shown to 
give reasonable results for several materials [12,49,50]. 

It is primarily in this context (creep/fatigue interaction) that the main 
distinctions between models I, II and III of Table II can be observed. For example 
utilizing the creep model of Table I (as will be the case henceforth) and the fatigue 
model I or III of Table II for the functions f and fp, respectively, one can 
immediately see that interaction between the two damage measures occurs from the 
outset. However if the function fp is comprised of Model II, interaction between 
creep and fatigue damage measures does not occur until crack initiation has occurred. 

4.0 CDM MODELS FOR INITIALLY ANISOTROPIC MATERIALS 

The Continuum Damage Mechanics approach, which supports the above 
described global measure of fatigue damage in monolithic isotropic materials, presents 
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two deficiencies which are inconsistent with a continuum approach. They are, i) the 
surface character of fatigue damage and ii) the fact that during the propagation stage 
the number of large defects is small. Interestingly enough, however, these deficiencies 
do not appear to have diminished the correct description of many experimental 
results in terms of cumulative fatigue damage. 

It is the authors contention that these two theoretical deficiencies will be 
removed when one establishes a representative volume element (RVE) on the 
mesostructural scale of a composite material; since now the surface character of 
fatigue will transcend to the interface region of each constituent (or phase) and the 
number of large defects occurring during the propagation period will similarly 
increase. This hypothesis is borne out by experimental observation [55,56] under 
specific thermal and mechanical loading histories. Although, there are a number of 
surface and environmental effects which will not necessarily transcend to the fibers, 
for example those effects which are associated with oxidation in the SiC/Ti 15-3 
system [57]. A nonrigorous definition of an RVE for a given composite material is 
given in Fig. 3 and in short insists upon the inclusion of a sufficient number of unit 
cells within the RVE to allow (statistically) the homogenization of the heterogeneous 
nature of the material. 

Figure 4 depicts, the extension of the previously described creep, fatigue and 
creep/fatigue damage mechanisms of Fig. 2 to the mesostructural scale, in which the 
surface or interface of a constituent (the fiber) plays the role (on the mesostructural 
scale) of a grain (and therefore grain boundary) on the microscopic scale. Here 
similar damage mechanisms are postulated to occur on the mesostructural scale as 
those on the microstructural scale, and due to the internal structure of the material, 
a sufficient number of defects will be present to allow a theoretically consistent 
continuum representation of creep and fatigue damage. 

In addition to the matrix cracking and formation of micro-cracks around the 
fibers, other damage modes must be included such as fiber breaking, fiber /matrix 
debonding and interlaminar cracks. It has been noted that damage occurs in the 
form of different multiple cracking modes and that there is no isolated single crack 
that dominates the development of damage. Similarly in composite laminates damage 
develops along preferred orientations, for example, matrix cracks in off axis plies 
typically are channeled by the fibers in those plies and interlaminar planar cracks 
grow along fibers in the neighboring plies [58]. 

These direction oriented damage modes suggest the need for some direction 
dependent damage variable (e.g., a vector or tensor representation). It is assumed in 
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the present work that the damage measure is a scalar but that the damage evolution 
is anisotropic. This is believed to be justified based on the strong initial anisotropy 
of the composite. Furthermore, as discussed earlier the definition of effective stress 
and the concept of remaining life will be utilized to measure and interpret damage. 

4.1 Transversely Isotropic Creep Damage 

Figure 4a is the mesostructural counter part to Fig. 2a and recently Robinson 
et al. [59] extended the creep damage model proposed by Leckie [40,41] and Leckie 
and Hayhurst [8,13] to metallic composites. This extension was accomplished by 
introducing into the isochronous damage function appropriate stress invariants that 
correspond to the local maximum transverse tension (stress normal to the local fiber 
direction) and longitudinal shear (shear stress on planes containing the fibers and in 
a direction along the fibers) within a unidirectional metallic composite. These 
invariants are included based on the anticipation that the associated stress may 
strongly influence void growth at the fiber-matrix interface (as this interface is 
postulated to play a role, on the mesostructural scale, analogous to that of grain 
boundaries on the microstructural scale and interfacial degradation); and, consequently 
may correlate with a creep rupture mechanism based on interfacial degradation 
through diffusion related void growth [56,60]. 

The physically meaningful invariants I-p ^ and Ig, representing the maximum 
transverse shear stress, longitudinal shear stress and the maximum normal stress in 
the fiber direction, respectively, are defined as (cf. [61,62]); 

1 

( 6 ) 


in which 
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J 2 “ 2 S ij S ij 


I = D.. S-. 
1J iJ 


1 =Vjk s ki 


D.. = d-d- 
iJ i J 


S* • U' • n 1 S' ' 

ij ij 3 kk ij 


and d- (i=l,2,3) are the components of a unit vector denoting the local fiber 
direction. An additional invariant is required for the present formulation (i.e., the 
maximum transverse tensile stress) and is given in [59] as: 


— < 2 + S > 


( 7 ) 


where : 

< f > 


= 0 if f < 0 
= f if f > 0 


A = 

A = D ij Of 


S = vTT 


l 


Thus the isochronous failure surface (A) is assumed [59] to be a linear function of 
the two invariants Jf and g 5? i.e., 


A = 


( + OL &) 


( 8 ) 


where : 
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in which a c is a material constant and a Q is a normalizing reference stress. An 
alternative isochronous failure surface (A), assumed here, can be written as; 

A = 77- &W + 3(l-V^c) S + “c 

in which the hydrostatic state of stress, maximum transverse shear and longitudinal 
shear stress are assumed to be the dominant damage measures. Note that the two 
isochronous surfaces defined in equations (8) and (9) are similar, provided /? c is taken 
to be zero, since S is related to the maximum transverse stress Jf (cf. eq.(7)), i.e., 

s = {*- \ (4 - 4)} 

The specific functional form for A must be determined experimentally as discussed by 
Robinson et. al. [59] and Leckie [40,41], 

The damage evolutionary law [8,40], extended now (through the isochronous 
failure surface) to account for initially transversely isotropic material symmetry, is 
taken as : 

dD c 

TJ- = C A r (1 - D c ) m (10) 

in which : 


n A 


m = 


A— 1 


- 1 


C = 


( A— 1 ) 
n A t r, 


and n, r, A and tg are material constants, as specified in [40,59]. For instance A is 
often referred to as the creep damage tolerance and measures the ability of the metal 
to withstand local straining at points of high strain concentration in a structure. 
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For further information regarding the characterization of this growth law the reader 
is referred to [59]. 

If one assumes the isotropic limit, i.e., D. . = 1/3 since D~ = 1; then, 


I = y Sjj = 0 


( 11 ) 


1 = 3 Vij 


and thus 


S = 'l = 6 S ij S ij “ 3 J 2 


i if ' 1 = I, 


J S ij S ij = 3 J 2 


!g - 0 


and 


1 1 

* = T a ii + 7T ^2 

thereby indicating that equation (8) does not reduce to the equation put forth by 
Leckie and Hayhurst, that is, 

A = 3^~ I A + ( 1+a c 1 ( 12 ) 

while equation (9) does reduce to the previous equation [8,40], if a = 0, such that, 

V 

A = -57 Vc 4 + wtfpi ) (is) 

as seen by comparing 3> in table I to the above. 
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4.2 Multiaxial Transversely Isotropic Fatigne Damage Model 

Figure 4b is the extension of Fig. 2b and illustrates the damage, i.e., the micro 
and mesocracks, one might encounter on the mesostructural scale when a composite 
material is subjected to fatigue type loading. The approach taken to model this 
damage will be the stress based approach discussed earlier. Extension of these 
previous models (Table II) is accomplished, once again by introducing anisotropic 
damage surfaces with appropriate invariants [59,61,62] that represent stress states that 
are likely to strongly influence the various damage modes in metallic composites. 

For instance, we assume that i) the transverse shear stress (1^) and thus implicitly 
the transverse normal stress, will dictate matrix cracking; ii) the longitudinal shear 
stress (I2) dictates interfacial degradation; and iii) the maximum normal stress (I,j) 
in the fiber direction will dictate fiber breakage. In this way, it is believed that the 
following fatigue limit and static fracture surface contain the most pertinent damage 
producing stress measures. Similarly, a normalized stress amplitude measure can be 
defined. That is, let 


( 4 j- 1 ) 

F 0 = J TTI {( 4 *'( l ) -m i +— ^ 1 2 + 9/4 I 3} (14) 

where the fatigue limit surface 2 is defined as ; 


1 

F f< - 1 = 2 m f - "ijM - 1 


(15) 


the static fracture surface is; 


1 - F u = 1 - max P (u) (ay(t)) 


(16) 


2 The authors acknowledge that not all monolithic materials exhibit a maked fatigue 0 
endurance limit, e.g., nonferrous alloys, however it is assumed (based on experimental 
evidence [63]) that a sufficiantly large class of composite materials will posses such a 
limit. Furthemore, we assert that fatigue lives beyound some specified value, say 
7 

10 , can be considered infinite for most practical design puposes and thereby allow 
the assumption of an endurance limit. 
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and the normalized stress amplitude is; 


F m = F m( <r ij) = 2 max [ F (m)( <T i/ t ) “ (17) 

wherein the (.) in eq.(14) is to be replaced by either £1, u, or m 
with 


U L - a U L • • 1 il - T uJ T iL i ' U L - ’ "u “ U l/ U T ’ 

\ = T u L / r u T ' m L = M L ’ u m = M l/ M T and = T m [ / T n> T 

being material parameters, some of them varying with the mean stress, that describe 
the initial anisotropy of the surfaces and m^ is defined by equation (40). Note that 
the subscripts L and T denote, respectively, the longitudinal (in the fiber) and 
transverse (normal to the fiber) directions. Now with these two surfaces and 
normalized stress amplitude it is a straight forward task to extend the models of 
Table II to account for materials which posses initial transversely isotropic material 
symmetry. The new multiaxial representations are now given in Table IV. 

One additional modification, besides the inclusion of a preferred direction, has 
been introduced into the fatigue models of Table IV and that is a hydrostatic stress 
measure (P) that combines both the Sines [64] and Crossland [65] criteria, i.e., 

P = ^ °Hmax + “ 0 °Havg. ( 18 ) 

This alternate hydrostatic measure (P) is felt to be required since a judgment as to 
the most applicable criteria is unavailable due to the lack of experimental data on 
metallic composites. In Table II, only the Sine’s criteria was considered, i.e., £ = 0. 
Finally, if one assumes that the ratios co^ tj^ a^, 77^, a> m , r) m , all are equal to one 
then the isotropic models given in Table II are easily recovered with only a slight 
modification to the coefficient ’a’ in models I and III and Cq in model II. 
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4.2.1 Uniaxial Simplification of Model I 


Due to the similarity in structure of models I, II and III and the fact that 
model I is the simplest to describe, only this model will be considered subsequently. 
Considering a uniaxial stress state (a) and a fiber direction oriented in the X-Y 
plane at an angle # from the X axis, i.e., dj = (cos#, sin#, 0). The following 
simplification of model I, cf. Table IV, is obtained; 


2 

a 

J 2 = ~Z~ 

I 2 = ^ (2 cos 2 # - sin 2 #) 2 a 2 

I = g- (4 cos 2 # + sin 2 #) a 2 

with 

*1 = * 2 A 1 

X 2 = A 2 

4 2 

*3 ~ 9 a A 3 
such that 


(19) 


( 20 ) 


dD„ = [H1-D./+ 1 ] 0 (Ccfl 2 


rVu-'iC 


m ' l( i-Dr.) 


dN 


with 


a = 1 - a 


<jcc u ( * max -<0-i> i 


<1—JCC a > 
v u max 


( 21 ) 


( 22 ) 
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where 


CC = 
m 


{( 4a m -1 ) A 1 + 


(4w 2 -l) 
v m } 


A n + Aq} 


'm 


(4«1 -—1) 


1 l 

cc u = 71 : A i + 


a- 


3 A 2 + A 3^ 


’ f / 


( 23 ) 


cc u = 


1 C 4 ^- 1 ) 

a 12 {( 4w u _1 ) A i + ^ A 2 + A 3^ 

U T U 


14 1 1 

Aj = -g - g- cos^0 - -g- sin^0 + g-g (2 cos^0 - sin^0)^ 


4 11 

A 2 = g - cos^d + g- sin"*# - g- (2 cos^0 - sin^0)^ 


(24) 


1 _ 2 /i\ 2 


Ag = (2 cos 0 - sin p) 


and henceforth, 


< f > 


= 0 if f < 0 
= f if f > 0 


with Mj j (P), ^P), a (P) and the oj and rj ratios are independent of (7g since the 

L 

modification factor cancels. Note once again that if u>q and t/q are equal to 1 (i.e., 
isotropy) then equations (21) and (22) become. 


dD p = [H1 -d/ + 1 ]“ lM L m (l-D F ) l' ? dN 


(25) 
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and 


(26) 


- ^ ff ma x a f l ^ a ^ / 
a = 1 - a 

/a -a ) 

\ u max / 

with a^a) = <r f ^(0)+(l - b a^O)) a 

and a = a 

which are identical in form to the uniaxial ONERA model I of Table III. 

Integrating equation (21) for N cycles, where «r and a are held fixed leads 
to : 



[1-(1-D F )^ +1 ] (1_o) 


1 + 0 ) (1-a. 


[VCC m (a- 3)F* 


If we cycle the material to failure (Np), which is defined to occur when Dp — 1, we 
find that 



Now rewriting equation (27), an expression describing the accumulation of 
damage, for a given loading cycle, in terms of the remaining life ratio can be 
obtained, that is : 

Dp = 1 - [1 - (N/Np) 1 /( 1_a V^ 1+ ^ ( 29 ) 

Note that the accumulation of damage is nonlinear and is also a function of fiber 
orientation, as is the number of cycles to failure. 
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Similarly if a two-level loading test is considered, one obtains, by integrating 
equation (21) in two steps, the same expression as in the case of isotropic materials 
i.e; 


N 2 /N p = 1 - (Nj/Np )P P = (1 - OjVfl - Qj) (30) 

2 1 

Except here o^, a Np and Np are all dependent upon the orientation of the 

1 2 

fibers (or conversely the load). Note that and N F are the failure lives (on the 

*1 * 2 

S-N curve) for the first and second loading conditions, respectively; while is the 
number of cycles applied at the first loading level and Ng is the remaining life at 
the second level. 


4.2.2 Nondimensional Uniaxial Parametric Study 


In order to conveniently examine the behavior of this fatigue model in detail 

for a variety of composite systems the following dimensionless form has been taken. 

Here both the maximum and mean stress have been normalized relative to the 

ultimate, static fracture, stress (a ) in the fiber direction, i.e., 

U L 


S 


max 


and 


max 



5 = 



(31) 


, respectively. Equations (21) and (22) can then be rewritten as: 


dDp 

dN 


[l-(l-Dp)^ +1 ] a (tftfj# 2 


.Ti Dp ) 


(32) 
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with 




a = 1 - a 


where 


[1-Vtftr s ] 

1 y u max 1 


m 


1 

JU A 1 + tT t A 2 + A 3^ 


(33) 


~ ~W fr 1 ) A i + 


( 4w f z“i) 


f 1 


An + Aq} 


(34) 


W u 


( 4 w*-l ) 

= {(4^-1) Aj + A 2 + A 3 } 


M, 


Jt 


it 


91 = 


(35) 


and Aj, and A^ are those defined previously. Again integrating equation (32) one 
obtains expressions relating, i) the number of cycles to failure relative to the applied 
normalized stress, 


[ \j #’#’ m ( S m a x "S") ] ^ / [ 1 V tftfu S m a x ] Wo^n 

wn — ( j 7w~ (s..-a,-| v ’ 

and ii) the damage accumulation relative to the ratio of remaining life: 


D p = 1 - [1 - (N/N F ) 1 /( 1_a )] 1 /( 1+/? ) 


(37) 
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4.2.2. 1 Direction Independent Parameters 


First, considering an isotropic material, i.e., Wq = t/q = 1, the effect of 

varying the parameters A A P> a, and the mean stress (5) on the S-N curve will 
be examined. To simplify and yet not limit the study, the parameters <r^ and 
will be taken to be independent of hydrostatic stress, that is to say b and b’ (given 
in Table IV) are zero. The baseline values assumed are those corresponding to 304 
stainless steel (see Table V and Fig. 5) and have been obtained from Lemaitre and 
Chaboche [12]. Due to the assumption of material isotropy, equation (36) can be 
simplified to: 


rS 


N F “ 


TT+7TT 


max 


Jl 


-ST 


P / [ 1 S max, 


max 


-(T+a) ] 


■> 


(38) 


Now as suggested by equation (38), if either the mean stress (?) or fatigue 
ratio ( are modified, the load level at which an endurance limit is reached is also 
changed. This is clearly shown in Figs. 6 and 7, where mean stress values of, 

? = 0.0, 0.2, 0.4 and 0.6 and fatigue ratio values of 91 = 0.1, 0.25, 0.4, 0.55 and 

0.7 are considered, respectively, while all other parameters are held fixed at the 
baseline values given in Table V. In both cases, as the value is increased so is the 
endurance level. 

Translation of the S-N curve is achieved, as indicated by equation (38), by 
varying the parameter Jl (i.e., Jt = 1.125, 2.25, 4.5 and 9.0) and ’a’ (i.e., a = 0.01, 
0.1, 0.5, 1.0, and 2.0) as shown in Figs. 8 and 9, respectively; while again holding all 

other parameters fixed. It is clear that increasing Jl shifts the S-N curve to the 

right (increases the number of cycles to failure at a given load) while increasing ’a’ 
shifts the S-N curve to the left (decreases the number of cycles to failure at a given 

load). Thus only the product (a Jl and not the individual values of a and A 
is important in determining the life at a given load level. However, the actual value 
of ’a’ does strongly affect the damage accumulation and therefore if any other type 
of damage (e.g. creep) is present , both parameters Jl and ’a’ need to be determined 
explicitly. 
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Finally, equation (38) indicates that modifying the parameter /3 will both 
translate the S-N curve and affect the abruptness of the transition from infinite life 
to immediate failure. Clearly, ft equal to zero is the lower limit and will give the 
most abrupt S-N curve for a fixed set of material parameters. Figure 10 illustrates 
the effect of varying /?, were /? = 0.5, 1.0, 3.0, 5.0 and 7.0 are considered. 

4. 2. 2. 2 Parameters Defining the Initial Transverse Isotropy. 

Considering the assumed baseline transversely isotropic material parameters 
given in Table V, the modification of the S-N curve with fiber orientation (or 
equivalently, load orientation relative to a fixed fiber orientation) is examined in Fig. 
11. As one might suspect, (as a composite material is designed to be stronger when 
loaded parallel to the fiber direction) both the static fracture strength and fatigue 
limit are reduced with an increase in fiber orientation. The amount of decrease in 
load carrying ability is clearly dependent upon the "degree" of anisotropy and is 
represented in this fatigue damage model by the u and tj ratios in the ultimate 
stress ( ## u )> fatigue limit ( tf#^) and normalizing stress amplitude ( && m )> see 
equations (36). Note that in the baseline material parameter set, all three 
anisotropic functions are considered to have the same degree of initial anisotropy. In 
this section an examination of the impact of perturbing the various anisotropic 
measures will be undertaken. 

Equations (34-36) and (24) indicate that the angle dependency is manifested 
through the functions Ap A 2 , and A^ which are associated with the longitudinal 
shear, transverse shear and normal stress components, respectively. This angle 
dependency is clearly shown in Fig. 12. At 0=0° both the longitudinal and 
transverse shear components are zero; at 0=45° the longitudinal shear stress 
component is a maximum and the transverse shear and longitudinal normal stress are 
equal; at 0=57° the longitudinal normal stress is zero and at 0=90° the transverse 
shear and longitudinal normal stress are equal. Figure 12 clearly indicates that when 
0=0°, changes in any of the measures of anisotropy (i.e., w’s or rf s) will have no 
affect since both the Aj and A ^ expressions are equal to zero. Similarly, when 
0=90°, any change in the shear "strength" measures (i.e., rf s) will not have any 
impact in the resulting S-N curves. Thus suggesting, that initial characterization of 
material parameters a, b, /?, # and JC should be conducted with longitudinally 
reinforced and loaded specimens. 
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4. 2. 2. 3 Variation of cj , (J r/ and u; . 

u’ it m 

Figure 13 illustrates the effect of varying the ultimate strength ratio 
(w =a I a ) on the S-N curve corresponding to a fiber orientation of 0=15°. The 

11 llj-^ 

values of u) are taken to be 2, 4, 5, 8 and 16, while all other parameters are held 
fixed at the baseline values given in Table V. As one might expect increasing the 
ultimate strength ratio affects only the low cycle fatigue portion of the curve (by 
decreasing the static fracture stress) while the endurance limit remains unaffected. 
Conversely, increasing the fatigue limit ratio (i.e., = 2, 4, 5, 8 and 16), while 

holding all other parameters fixed, decreases the endurance limit while leaving the 
ultimate strength unchanged. See Fig. 14 for an example at a fiber orientation of 15 
degrees. The impact of changing from 5.0 to 3.0 for fiber orientations of 0, 15, 
30, 60 and 90 is shown in Fig. 15. Comparing Figs. 11 and 15 one observes a 
similar trend (as that observed in Fig. 13) for all angles except 0=0° , which remains 
unaffected for reasons indicated earlier. Note that although the present model 
provides significant flexibility, in that and can be varied independently; in 
reality one would expect a relationship to exist between the ultimate strength and 
the endurance limit, thus placing restrictions on the ranges of these parameters. 

Furthermore, by varying the normalizing stress ratio o> m , a horizontal 
translation of the transition portion of the S-N curve can be obtain while the 
ultimate stress and fatigue limit remain unchanged. An example is shown in Fig. 16 
where takes on the values of 2, 4, 5, 8 and 16, assuming a fiber orientation of 
15 degrees, while all other parameters once again remain fixed. Clearly, increasing 
the ratio shifts the S-N curve to the left (decreasing the number of cycles to failure 
for a given load). 

Thus, by merely assuming different degrees of anisotropy for the three functions 
if if , ifif-4 and if# , i) either end of the S-N curve can be modified while leaving 
the other end unchanged, or ii) the center portion can be horizontally translated, 
while the ultimate and fatigue limit stresses remain unchanged. This suggests that 
the model has sufficient flexibility to fit a broad class of materials. 

The impact of varying the degree of anisotropy, for example on the S-N 

curve can be seen for all angles of fiber orientation by plotting the functions if iff ^ 

and V if if u , <J if if m _ In Fig. 17 only <J if if has been shown versus angle of 
orientation 0, as the other two functions will have similar behavior, except for being 
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scaled by a multiplying factor. Examining Fig. 17 it is evident that for a specified 


ratio the function V increases as one increases the angle 6 ; however, 
depending upon the magnitude of an angle can be found where this increase is at 
a decreasing rate. This explains why, for larger angles, one sees an incrementally 
smaller change in the response curve with increasing angle (i.e. , variation in the S-N 
curve), see Figs. 11 and 15. Finally, one more important point can be discerned 
from Fig. 17, and that is the need for the ratio u to always be greater than or 
equal to 77 since if u> < 77 a minimum occurs in the function at an angle other than 
zero, which violates physical reasoning. 

4. 2. 2 . 4 Variation of 77 ^, 77 ^ and 77 m - 

Variation of the shear ratios 77^, 77^, and 77 m has a similar impact on the 
overall trends, but not actual magnitudes, of the S-N curve, as did their respective 
normal stress ratios. Examples involving the variations of each ratio (for values of 

1, 2, 4, and 5) with a fiber orientation of 15 degrees are shown in Figs. 18, 19 and 
20, respectively. The interaction between the shear strength and the normal strength 
ratio are shown in Fig. 21, for an a> u =5.0, and in Fig. 22, for an w u =16.0, when the 
77 u ratio takes on values of 1, 1.5, 2, 3, 4, and 5 in the case of Fig. 21 and 1.0, 1.5, 

2, 4, 8, and 16 in the case of Fig. 22. Note how in both cases (i.e., 07^=5. 0 and 
k> u =16.0) the shape of the function is greatly affected, particularly for low 
values of 77. Also it is clear that when 77=1 a maximum at some intermediate angle 
is obtained. The ramification of this intermediate maximum is unknown at this time. 
As with the cases involving the normal stress ratios, and , one would expect 
that a relationship would exist between 77 u and 77^ such that the practical ranges of 
these parameters are not independent. 

4.2.3 Required Exploration and Characterization Experiments. 

Here, the determination of the various material parameters will be discussed. 

Fatigue of Isotropic Material 

Considering first the case of isotropic materials, it is assumed that S-N curves 
(shown schematically in Fig. 23) exist for different stress ratios (R=a m in/cr ni ax) and 
that the ultimate static fracture strength (a u ) is known. Given this data, one may 
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then obtain a plot of fatigue limit (tr^) versus mean stress (cf), as shown in Fig. 24, 
from which the material parameter b can be easily found. Note that a linear 
dependence is assumed in Fig. 24; this, however, may not be sufficient to describe 
the material and one might have to consider an alternative representation. 

Rewriting the isotropic simplification of equations (23) and (28) and taking the 
natural logarithm of both sides, one can obtain the following expression : 


Kcr m o 

— 
o u — <j\ 


& 1 

N F 
m / r 


= In 


a(l +0) M 


- (9 lnfff.-F) 


(39) 


or 


Y„-/0X 


Now plotting Y versus X (see Fig. 25) and fitting a straight line through the 
data, the material parameter (3 is obtained from the determined slope and the 

product (a M - ^) may be determined from the intercept. Similarly solving equations 
(23) and (28) for M, i.e., 


M = [ n f (i-a)] 1 ^ 


M = (cr ra -a) [Np (1+/3) 


/( o m — a )— a f u 
\ & u (7 m / 


«Tf 1 J 


(40) 


It becomes apparent that M represents a normalizing factor such that it ensures a 

certain life (Np) for a specified stress amplitude, <T 3L = ( a m ~^)y see ^ig. 26. Clearly, 
M is also dependent upon the slope /? and the material constant a, thus leading one 
to the conclusion that the remaining life measurements are insufficient to completely 
characterize the damage evolution equation, i.e., the value of ’a’. Only the product 

(a M P) may be determined directly, by either, using the previous Y intercept (Yq) 
of Fig. 25, i.e.; 



1 

(1+jS) e Y ° 


( 41 ) 


or selecting a point on the S-N curve (see Fig. 23 or 26) such that 


25 



( 42 ) 


< a — a> a tt 
_0 a u m ft 

a M = N F ( 1 + ^) <ff a“‘ 7 f l> 17 u 

This indeterminacy (lack of independence between ’a’ and M) has been shown 
not to be important, as long as only fatigue damage accumulation is considered, even 
under complex loading histories. However, in the case, for example, of creep-fatigue 
interaction the complete characterization (independent measurement of parameter ’a’) 
of the model is required. 

To accomplish this, an indirect measurement of damage accumulation can be 
obtained from the change in the stress-strain response during the fatigue process 
using the effective stress concept. One method of experimentally measuring the 
damage accumulation is to monitor the change in elastic response (Young’s Modulus) 
with applied cycles of stress, that is 


E 

D = 1 - g (43) 

where E is the effective (current) modulus and E is the initial one. A number of 
researchers have used this procedure for both metallic and composite materials 
[15,31,34] and in fact such measurements have been observed to correlate well with 
quantitative micro-crack evaluations as well [42,43]. Now given D as a function of 
remaining cycles (N/Np), all that is required to determine ’a’ is to plot the natural 
logarithm of D versus that of N/Np, as shown in Fig. 27, calculate the average slope 
Z for various stress amplitudes, and solve for ’a’, i.e.; 


a 


< 1 + 0 ) <%-%> 

Z <a a — a f 


“it 


(44) 


With this, all material parameters, i.e., a, /?, b, a , a ^ and M, have been 
determined in the case of isotropic materials. 

If model II is selected, the function describing the influence of frequency needs 
to be characterized, that is: 
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This is accomplished by plotting semi-logarithmically N. versus (cr m - <7 ^)/( <7 u -lT m ) 
for a variety of frequencies; where N. is obtained by making use of the empirical 
expression put forth by Manson [66], that is, 


N. = N p - 14 N f 0 6 

Thus the intercept gives C(^) while the slope is Now plotting C{u) versus v 

and rewriting the expression for C(V), for example 


CM 



1 + 

V 0 ' 


V 



jVLy) 

O j/7 


and plotting C{v) versus C(u)/i/^ we see that and C 0 are the slope and intercept, 
respectively, given a value for 7. The point here is to plot a number of curves 
corresponding to different 7’s and selecting the one that best fits the data. This 
procedure is described in greater detail in reference [25]. 

Fatigue for Transversely Isotropic Materials 
For transversely isotropic materials, e.g., hexagonally packed unidirectional 
composite materials, the determination of the material properties follow a similar 
process, except now additional testing is required both in tension and torsion in the 
longitudinal and transverse directions, to define the extent of initial anisotropy, i.e., 
the ratios 77 , u ^ 77^, and w , 7? m . Figure 28 illustrates schematically the 
expected S-N curves for a composite material loaded longitudinally (0=0°) and trans- 
versely (0=90°) at a given mean-stress. Given such data, as well as the associated 
shear S-N curves, the various u and 77 ratios can be determined. Note that in the 
present theory /?, a and b are assumed to be independent of fiber orientation (i.e., 
dj). This assumption may need to be modified if experimental evidence suggests 
otherwise, therefore during the initial exploration and characterization of a class of 
composite systems the procedures illustrated by Figs. 24, 25, and 27 should be 
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repeated four times, that is for the case of normal and shear type loadings in the 
longitudinal as well as transverse directions. 

Clearly the above discussion suggests a significant, and potentially quite 

expensive, experimental program. This experimental program may however be able 

to be reduced and augmented through numerical simulation using a suitable 

homogenization technique with periodic boundary conditions. Furthermore utilization 

of the homogenization technique should allow micro-mechanical effects (e.g., bond 

strength, volume fraction, etc...) to be included in this phenomenological theory, 

through the parameters, , u> u , r? u , and 77^ 

L L 

One important micromechanical aspect which has been intentionally neglected in 
this initial formulation is the role of residual stresses (which are expected to be 
significant) in the life of metallic composites. However, at this time insufficient 
experimental data exist to verify the various extensions of the presently described 
models, let alone admit the introduction of further complicating factors. Therefore 
this aspect will be addressed in future work. 


Creep During Cyclic Loading for Isotropic Materials 


A final remark regarding the characterization of the required material 
parameters. The parameter r in the rate form of the delay stress equation in Table 
I, can be obtained most easily by assuming a sinusoidal loading history. For 
example in uniaxial form the required equation set is 

a = A sin wt 
X = C sin (iiA + <p) 

<J2L= - x )/ r 

cTF 

u = 2x v 

and the solution is 


X = 


(1 + tut) 


Tg — j sin ciA — tlj cos 


An expression relating the maximum delayed stress (X ) to the maximum applied 

max 7 


stress (a ), 
v max'’ 
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G 

max 

■^max ~ . . 22^0.5 

( 1 + u t ) 

and is obtained by finding the stationary point of the above equation. Now 

assuming that X = <r at a sufficiently high frequency of loading , an 
max mean 

expression defining r is obtained 

r = — ^ — { (a la — 1.0 } 

27T v 1 v max' mean' J 

Clearly, the suitability of this frequency dependence should be examined when the 
load history is a fully-reversed one. 

4.3 CREEP-FATIGUE INTERACTION 

Figure 4c is the extension of Fig. 2c and illustrates schematically the 
interaction of intergranular and transgranular defects. As for isotropic models it is 
surmised that the presence of cavities allow for easier crack propagation and that the 
increase in stress intensity at a crack tip causes an increase in the nucleation and 
coalescence of voids. 

Analytically this interaction is again represented using the effective stress 
concept, by assuming that the mechanical effects of creep and fatigue damage can be 
directly added, (see equation (5)). Thus the two functions (f and fp) can be 
determined independently from pure tensile creep tests and pure high frequency 
fatigue tests. 

It is primarily in this context (creep/fatigue interaction) that the main 
distinctions between models I, II and III of Table IV can be observed. 

5.0 CONCLUSIONS 

Three isothermal multiaxial continuum damage mechanics models for the creep, 
fatigue, and creep/fatigue interaction of a unidirectional metal matrix composite 
volume element have been presented. The intended applications are reinforced 
structures in which the fiber direction may vary throughout but a single fiber 
direction can be identified locally (local transverse isotropy) within a given volume 
element. Each model is phenomenological, with varying degrees of complexity to 
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accurately predict the initiation and propagation of intergranular and transgranular 
defects over a wide range of loading conditions. 

The development of these models are founded on the definition of an initially 
transversely isotropic fatigue limit surface, static fracture surface, normalized stress 
amplitude function and isochronous creep damage failure surface, from which both 
fatigue and creep damage evolutionary laws can be obtained. The anisotropy of each 
model is defined through physically meaningful invariants reflecting the local stress 
and material orientation. All three transversely isotropic models have been shown, 
when taken to their isotropic limit, to directly simplify to previously developed and 
validated creep and fatigue continuum damage theories. 

Results of a nondimensional parametric study illustrates i) the flexibility of the 
present formulation when attempting to characterize a large class of composite 
materials and ii) its ability to predict anticipated qualitative trends in the fatigue 
behavior of unidirectional metal matrix composites. Additionally, the potential for 
the inclusion of various micromechanical effects (e.g. bond strength, volume fraction, 
etc.), into the phenomenological anisotropic parameters have been noted, as well as a 
detailed discussion regarding the necessary exploratory and characterization 
experiments needed to utilize the featured damage theories. 

Two potential drawbacks to the present formulation are i) the scalar damage 
measure employed and ii) the expensive experimental program required. Future work 
in this area will include the examination of the utility of the present formulation to 
predict high temperature applications, the initiation of an experimental 
characterization and validation program and the numerical implementation into a post 
processing life prediction computer code. 
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APPENDIX A - NOMENCLATURE 


Stresses: 


a ij 

is 

the 

d 

a. . 

is 

the 

ij 

a 

a 

is 

the 

S-- 

ij 

S 

is 

the 

is 

the 

max 

3 

is 

the 

a 

is 

the 

a 

is 

the 

m 

ff a 

is 

the 

^Havg 

is 

the 

^Hmax 

is 

the 

a 

is 

the 

a o 

is 

the 


Cauchy stress tensor 

delayed stress tensor 

uniaxial delayed stress 
deviatoric stress tensor 

normalized uniaxial maximum stress 

normalized uniaxial mean stress 

uniaxial mean stress 

uniaxial maximum applied stress 

applied stress amplitude, i.e., a m -<7 
uniaxial mean hydrostatic stress 
uniaxial maximum hydrostatic stress 

uniaxial effective stress 
reference stress 


Invariants: 


J 


o 


Jf 

S 

A 

'A 


is the stress invariant representing the maximum principal stress 

is the first invariant of the total stress (hydrostatic stress) 

is the second invariant of the deviatoric stress (expressed for shear) 

is the second invariant of the deviatoric stress (expressed for normal 
stress) 

is the invariant representing transverse shear stress 
is the invariant representing longitudinal shear stress 

is the invariant representing the maximum normal stress in the fiber 
direction 

is the invariant representing maximum transverse tensile stress 
is the square root of 1^ 

is the square root of I„ 

is the invariant representing the total stress in the fiber direction 
represents the maximum amplitude in stress 
represents the avaraged mean stress 
represents the maximum octahedral stress 
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Material Parameters: 


a /? 
c, c 



c 

r,k 
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JC 

91 1 
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T il 
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oV a ol 
oi ,a oi 
op ,<T op 
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u 
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a, a 


b,b’,b’ 

a 

P 

e 

V 

u 


Creep Damage 

are the coefficients indicating the strength of contribution for the J Q and 

^ invariants, respectively. 

is the creep damage normalizing factor 

are the exponents in the creep damage evolution equation 
is related to the response time of the material 

Fatigue Damage 

is the ratio of stress amplitude constant to static fracture stress 
is the ratio of endurance limit to static fracture stress 
is the normal stress endurance limit 

is the ultimate normal stress or static fracture stress 

is the normalizing stress amplitude 

is the ultimate shear stress or static fracture stress 

is the shear stress endurance limit 

is the uniaxial initiation endurance limit 

is the uniaxial propagation endurance limit 

are the initial multiaxial and uniaxial endurance limits, respectively 
are the initial multiaxial and uniaxial initaition limit stresses, respectively 
are the initial multiaxial and uniaxial propagation limits, respectively 
is the ratio of longitudinal to transverse normalizing normal stress 
amplitude 

is the ratio of longitudinal to transverse ultimate normal stress 

is the ratio of longitudinal to transverse normal fatigue or endurance 

limit stress 

is the ratio of longitudinal to transverse shear static fracture stress 
is the ratio of longitudinal to transverse shear fatigue limit stress 
is the ratio of longitudinal to transverse shear normalizing stress 
amplitude 

is the scaling factor for the stress dependency in the fatigue damage 
exponent 

are parameters indicating the effect of mean stress 

is the fatigue damage variable exponent which can be a function of stress 

is the exponent on the normalized stress amplitude 

is the exponent on the normalized stress amplitude during initiation 

is the angle between the fibers and coordinate axis 
is the frequency of the applied loading cycles 
is the reference load frequency 

is the exponent of the frequency term 
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C Q is a pre-multiplying parameter in the initiation section of model II 
E Undamaged Young’s Modulus 

E Damaged Young’s Modulus 

( subscript denotes longitudinal properties 

( )rp subscript denotes traverse properties 


Miscellaneous: 


d i 

D-= 

>J 

*ij 

A 

D 

* 

D 

D c 

d f 


is the unit vector denoting the local fiber direction 


d i d j 


is the second order direction tensor 


is the Konecker delta function 


is the isochronous creep damage failure surface 
is the scalar variable representing damage 

is an alternate damage variable that can be used interchangably with D 
is the creep damage measure 

is the fatigue damage measure 


F^ is the fatigue limit surface 
F u is the static fracture surface 


F m is the normalized stress amplitude 
P is the combined hydrostatic mean stress measure 

( is the parameter in P, weighting the influence of the maximum 
hydrostatic stress 

N is the variable denoting number of cycles 

Np is the number of cycles to failure 

Nj is the number of cycles to crack initiation 

Np is the number of cycles for crack propagation 

u is the frequency of loading 

t is time 

T is temperature 

< > are the Maclaurin brackets 
H( ) is the Heavyside step function 

3> is the isochronous creep damage failure surface perposed by Leckie and 
Hayhurst 

A is the area on a face of a representative volume element 

A is the effective area on a face of a representative volume element 

Ap is the damaged area on the face of a representative volume element 

p is the phase shift in the sinusoidal load history 
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Table I : The Leclde-Hayhurst Creep Damage Model 


dD 


c 


dt 


Multiaxial: 


r x(<r ij)i 




with 

= “c A^ij) + J 2<' 7 ij) 

J 0 (° r ij) = m “ a u ; A^ ij) = "kk ; = 



S 


i j 


where 


y kk 

S- • = a- ■ n — S- ■ 

ij ij 3 ij 


, d 

d a ■ 


i J 


a ■ • —a ■ ■ 
i J i J 


d t 


and cr- j is replaced by cr^ . in the creep model 


A c , r, k, a c , fi c , and r are material parameters, A £ and r are typically taken to be 
temperature dependent. 


Uniaxial 


dD 


c 


dt 
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Table II : Fatigue Damage Models (ONERA NLCDM Multiaxial Form) 
Model I : 



^ S ^ ^ ^ 

where a = 1 - a 3 - 

<a u -^ 3 > 

M(^) = M q (1-b 

= (1— b’ 


Model II : 


Micro initiation : 




1 + 

V o' 


V 


< — S j ( ^2 ) 




< £r u -^ 3 > 


sM) = s oi f 1 - 1 *" 


Micro propagation : 



a = 1 - a 

y^> = y ^ ^ 
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Table II (Conclusion) 


Model m : 



^ S ^ ( 0^2 ) ^ 

where a = 1 - a — >— — 

< a u -^ 3 > 

M(^) = M Q (1 - bc^) 

S/^) = (1 - b’ 

S p(^) = S 0p (!- b ’^) 

Note M q> S^, S oi , S 0p , a, b ,b> ,b" , 0 , 7 and 
with a also taken to be temperature dependent. 

For all three models : 

1 

^ = 2 max [ max ( J 2 [^ i j( t i)- ff ij( t 2)])] 
l l l 2 

1 

^ 2=2 { max 4( a ij( t ))+ min ^(Oij^))} 

t t 

^ — max J 2 (^jj(t)) 
t 

where and Jg are defined in table I. 


] fi H(^-S p (^)) dN 


a are material parameters 


maximum amplitude 


mean stress 


maximum octhaedral stress 
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Table III : ONERA Fatigue Damage Models in Uniaxial Form 
Model I : 

d°F = IHl-D F )^ +1 ]° ImXT) ( X-Dp)!^ dN 
<c 

a = 1 - a — 

aft) = ff/0)+(l-b aft)) a 
M(<t) = Mg (1-b a) 


Model II : 


Micro initiation : 





< v. 


m 


< a 


~ a n 

U~V 


(a)> 

> 





Micro propagation : 

(T — (7 < o-—v , „ ( ~a 

O i i * in n HI t p 

dD p - [H 1_D p) ] [m ( a ) ( 1-D p )1 <a n~ a m > 

a = 1-a 

M(cr) = Mq (1-b cr)) 

a tp = < V 0 ) +(1_b V 0)) * 

where N P = N-+N 

r 1 p 
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Table m (Conclusion) 


Model m : 


dDp = [H1-D/ +1 ]“ ln (>) (i-P.)] " H ('nT‘ 


o—a 

m 


8 


<a m~ a t ( ) > 


a = 1 - a 


< a— a> 
u m 


= */°)+( 1 “ b a £°)) a 

M(o) = M q (1-b a) 


% = V° )+(1_b %(°)) a 

= 1 if f > 0 
= 0 if f < 0 


H(f) 


where, in the three models M Q , a^O), <j^( 0), <7^(0), b, a, 0 , 7, <r 
as material properties. 
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-*(*» dN 


are defined 


■ Table IV : Multiaxial Transversely Isotropic Fatigue Damage Models 
Model I : 


where 


and 


Model n 


and 


with 


dDp = [l-(l-Dj/ +1 ] a [ (l _p >] " dN 


m 




<F t 


a = 1 - a 


< 1 _F u> 


a U = a t l (°) ( 1_3b ’ p ) 

L L 


M l = M 0 (1— 3b P) 
L 


Macro initiation 
C 


0 


N. 


1 + 

V o' 

"7 

<i 


h V 






iT 


\i = Xi (0) (Mb " p) 


Macro propagation : 

*f< L P p 


<% — (Hi- D / +1 !“ [ (1 “ D j f — v - dN 




P' < 1 F U > 


a = 1-a 


\ P = Xp ( 0) (1_3b ’ p) 


M l = M 0 (1— 3b P) 
L 
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Table IV (Conclusion) 


Model in : 


dD F = "ttjp [ ( 1 — D p ) 1^ dN 


<F f<“ 1> 


with 


a = 1 - a 


<1 -F U > 


X = x (0) (1 " 3b ’ p) 
Xp = Xp ( 0) ' 1_3b ' P * 


M l = M 0 (1— 3b P) 
L 


where a, b, b\ b," p, 7 , <r f ^( 0 ), <7^(0), ^ L p(°)> M 0 L are material 

coefficients with a being taken to depend on temperature. 

U L 

Note that : 

1 

°Havg = T { max (° r ii( t ))+ m i n (^(O)} mean hydrostatic stress 
t t 


■'Hmax = XW) 

P = ^ <T Hmax + ( 1 ~^ ff Havg 
t - is time during a cycle 


max hydrostatic stress 

(combines the Sines and Crossland 
criteria) 
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Table V : Baseline Dimensionless Parameters for Isotropic and Anisotropic Materials 


Isotropic 


Transversely Isotropic 

St 

— 

0.25 

SI 

— 

0.25 

Jt 

— 

4.5 

Jt 


4.5 

0 

= 

3.0 

0 

= 

3.0 

a 


0.1 

a 

= 

0.1 

5 

— 

0.0 

3 


0.0 

w 


1.0 

U) 

= 

5.0 

m 



m 



“U 

= 

1.0 

U u 

=z 

5.0 

u„ 

— 

1.0 

U) 

= 

5.0 

u 



u 




= 

1.0 


= 

2.0 


= 

1.0 


= 

2.0 

% 

= 

1.0 

\ 

= 

2.0 



Volume element 



Virgin 



Damaged 


Equivalent 

virgin 


Figure 1 . — Schematic of effective stress concept and 
equivalence in strain. 
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crack 

| 

propagation | 

crack i 

propagation 
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i 

initiation 
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-L 


initiation ■ 



0.01 

Classical 

crack 

initiation 


0.1 1 

Present definition 
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mm 


Damage mechanics 


1 „ Fracture 
I mechanics 



(a) Creep damage: coalescence of cavities and (b) Fatigue damage: nucleation of slip bands, microcracks and trans- 

intergranular defects. granular defects. 



(c) Creep/fatigue damage: Interaction of intergranular and 
transgranuiar defects. 


Figure 2. — Schematic of different damage modes and the associated scale In a monolithic metal. 
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Figure 3. — Nonrigorous definition of representative volume 
element. 


Micro crack . Mesostructural 
propagation crack 

' Initiation ' 

I I 


Mesostructural . 
crack j 

propagation * 


Macro . Macro crack 

crack J propagation 

initiation * 
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0.1 


mm 



Meso 

initiations 



(a) Creep damage by way of Intergranular 
void growth. 


(b) Fatigue damage: initiation and propagation of 
transgranular defects. 


Stress 

Meso- concen- 



(c) Creep-fatigue: interaction of Intergranular 
and transgranular defects. 


Figure 4.— Schematic of different damage modes and the associated scale in a metallic composite. Fibers on 
mesostructural scale are analogous to grain boundaries on microstructural scale. 
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Figure 5.— Illustration of an S-N curve for the dimensionless 
isotropic baseline material parameters of Table V. 
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Figure 6.— Illustration of the effect on the S-N curve of vary- 
ing the mean stress S = 0, 0.2, 0.4, 0.6, while holding all 
other parameters fixed. 






Figure 7.— Illustration of the effect on the S-N curve of vary- 
ing the fatigue ratio = 0.1 , 0.25, 0.4, 0.55, 0.7, while hold- 
ing all other parameters fixed. 


Figure 8. — Illustration of the effect on the S-N curve of vary- 
ing the parameter M - 1 .1 25, 2.25, 4.5, 9.0, while holding 
all other parameters fixed. 




Figure 9.— Illustration of the effect on the S-N curve of vary- 
ing the parameter a = 0.01 , 0.1 , 0.5, 1 .0 and 2.0, while hold- 
ing all other parameters fixed. 


Figure 10. — illustration of the effect on the S-N curve of vary- 
ing the parameter 0 * 0.5, 1 .0, 3.0, 5.0 and 7.0, while hold- 
ing all other parameters fixed. 
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Figure 11 .—Illustration of the degradation in load-life re- 
sponse with variation In fiber orientation, 0 = 0, 1 5, 30, 60, 
and 90 degrees. 


Figure 1 2. — The angle dependency of the longitudinal shear 
(A^, transverse shear (A 2 ) and normal stress (A 3 ) compo- 
nents. 




Figure 13. — Illustration of the effect on the S-N curve of vary- 
ing the ratio u) u = 2, 4, 5, 8, and 16, at an angle of 1 5 de- 
grees, while holding all other parameters fixed at the base- 
line values. 


Figure 14. — Illustration of the effect on the S-N curve of vary- 
ing the ratio = 2, 4, 5, 8, and 16, at an angle of 1 5 de- 
grees, while holding all other parameters fixed at the base- 
line values. 



Figure 1 5. — illustration of the effect on the S-N curve of vary- 
ing the angle of fiber orientation 0 » 0, 1 5, 30, 60, and 90 
degrees with the ratio tu u = 3 and ail other parameters fixed 
at the baseline values. 



Figure 16. — Illustration of the effect on the S-N curve of vary- 
ing the ratio u> m = 2, 4, 5, 8, and 16, at an angle of 1 5 de- 
grees, while holding ail other parameters fixed at the base- 
line values. 
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Figure 1 7. — Shows versus angle of orientation 0, for vary- 

ing ratios of u> u , i.e., 1.0, 2.0, 3.0, 5.0, 8.0, and 1 6.0, when = 
2.0 and all other parameters are those in table V. 



Figure 19. — Illustration of the effect on the S-N curve of vary- 
ing the ratio = 1 , 2, 4, and 5, at an angle of 1 5 degrees, 
while holding all other parameters fixed at the baseline 
values. 



Figure 1 8. — Illustration of the effect on the S-N curve of vary- 
ing the ratio ti u = 1 , 2, 4, and 5, at an angle of 1 5 degrees, 
while holding all other parameters fixed at the baseline 
values. 



Figure 20. — Illustration of the effect on the S-N curve of vary- 
ing the ratio t| m = 1 , 2, 4, and 5, at an angle of 1 5 degrees, 
while holding all other parameters fixed at the baseline 
values. 




Figure 21 . — Variation of <€<€ u with angle of orientation 0, for vary- Figure 22.— Variation of with angle of orientation 0, for vary- 
ing ratios of tj u , i.e., 1 .0, 1 .5, 2.0, 3.0, 4.0, and 5.0, when o» u = 5.0 ing ratios of t] u> i.e., 1 .0, 1 .5, 2.0, 4.0, 8.0, and 1 6.0, when to u = 
and all other parameters are those in table V. 1 6.0 and all other parameters are those in table V. 
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Figure 23. — Schematic of typical S-N curve for an Isotropic 
metallic material, assumed to be known experimentally. 
Assuming R * ^ m jr/ CT max anc * the mean stress <r to be held 
constant. 



Figure 24. — A plot of fatigue limit (oy) versus mean 
stress (5), from which ct^{ 0) and b can be obtained. 




Figure 25. — A log-log plot of scaled cycles to failure versus Figure 26. — A schematic of an S-N curve defining the mean- 

stress amplitude is shown, from which the material parameter Ing of the normalizing factor M. 

p can be found. 



Figure 27.— A log-log plot of damage (D) versus remain- 
ing life (N/N F ) for various stress amplitudes. 



iog (N) 


Figure 28. — A schematic depicting the expected S-N curves 
for a unidirectional metallic composite, when loaded longi- 
tudinally (0 = 0.0) and transversely (0 = 90). 
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Additionally, the potential for the inclusion of various micromechanical effects (e.g. fiber/matrix bond strength, fiber 
volume fraction, etc.), into the phenomenological anisotropic parameters is noted, as well as a detailed discussion 
regarding the necessary exploratory and characterization experiments needed to utilize the featured damage theories. 
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